Geometric phases for AMevel systems through unitary integration 
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Geometric phases are important in quantum physics and are now central to fault tolerant quantum 
computation. For spin-1/2, the Bloch sphere S 2 , together with a U(l) phase, provides a complete 
SU(2) description. We generalize to JV-level systems and SU(iV) in terms of a 2(N — l)-dimensional 
base space and reduction to a (N — l)-level problem, paralleling closely the two-dimensional case. 
This iteratively solves the time evolution of an A^-level system and gives (N — 1) geometric phases 
explicitly. A complete analytical construction of a S 4 Bloch-like sphere for two qubits is given for 
the Spin(5) or S0(5) subgroup of SU(4). 
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Coupled quantal systems with N states, and their 
time evolution are of interest in broad areas of physics. 
Laser coupling between three or more atomic or molecu- 
lar states, oscillations between the three flavors of neutri- 
nos, and logic gates in quantum computation (or cryp- 
tography and teleportation) are but a few examples of 
the widespread interest in few-level systems. In a se- 
ries of papers [1] with specific applications to three- and 
four-level problems of quantum optics and quantum in- 
formation, we have explored a semi-analytic technique 
of "unitary integration" [2] for the time dependent op- 
erator equations involved. For Hermitian Hamiltonians, 
unitarity is preserved at each step, no matrix inversions 
are involved, and often the problem reduces to solving a 
single or a small set of Riccati equations for classical func- 
tions. The method also extends readily to non-unitary 
evolution with dissipation and decoherence [1, 3] that 
occurs in quantum computation. The evolution operator 
is written as a product of exponentials, each exponent 
containing one of the operators of the algebra and a mul- 
tiplicative classical function of time, these functions cho- 
sen appropriately to solve the time-dependent equations 
of interest. 

The increasing number of functions and coupled equa- 
tions that they obey (for m qubits, N = 2 m ) may 
have discouraged more applications of unitary integra- 
tion. This Letter presents a solution by giving a system- 
atic and compact scheme for general N, in a stepwise 
reduction that parallels the N = 2 single spin case. We 
also connect to another area of research, the study of 
phases that are central to the fields of quantum informa- 
tion and quantum computing. "Geometric" phases that 
depend only on global features [4] may be especially im- 
portant for quantum computing, possibly providing high 
fidelity and fault-tolerant operation [5]. Berry's initial 
discussion in terms of an adiabatic evolution and a single 
phase described by U(l) symmetry have been generalized 
extensively [6]. Mathematicians describe it as a holon- 
omy on a fiber bundle [7] . 

What has been missing, however, is a simple descrip- 
tion of the parameters and, specifically, phases for a N- 
level problem of the sort available for two levels. For 
the group SU(2) of a single spin, or qubit, two coordi- 



nates of a "Bloch" sphere S 2 [8], together with a U(l) 
phase, provide a complete description. In the language 
of differential geometry [9], SU(2) is thereby viewed as a 
fiber bundle of the base manifold S 2 and the U(l) fiber 
[7]. Generalizations to larger groups, especially SU(N) 
groups that describe higher spins or multiple qubits, are 
desirable but are less straightforward [10]. For instance, 
a two-qubit system for logic gates [11] has SU(4) sym- 
metry. Bloch vectors for N = 3,4, and general N [12] 
have been given but not the phase aspects and explicit 
evaluation of the evolution operator that we now provide. 
While having much in common with more mathematical 
treatments in [10, 13], our emphasis is on concrete con- 
structions through nothing more complicated than ma- 
trix algebra. 

The phase of a state's wave function being only ac- 
cessible in comparison with another reference state, we 
find more useful the "operator- valued phase" , which may 
be called the "W-phase", of the evolution operator U(t) 
since it is referenced to the unit operator at some ini- 
tial time. The continuous connection to the starting unit 
operator also means that there are no 2n7r ambiguities 
in these "non- modular" phases [10, 14]. Our method of 
"unitary integration" [1, 2] expresses U as a product of 
exponentials, each with one of the (N 2 — 1) generators in 
the exponent. Thus, for SU(2), U is a product of three 
exponentials in the Pauli matrices, cr_|_, <7_, cr z , the last of 
which provides the W-phase which can be further divided 
into dynamical and geometrical parts. For SU(iV), we 
construct a similar product form of two factors involving 
nilpotent exponents for the base and a third factor that 
is block diagonal in (N ~ n) x (N — n) and n x n blocks 
for the fiber. 

The use of nilpotent operators is central to our deriva- 
tion. The complex time-dependent coefficients multi- 
plying each operator in an exponent obey coupled first- 
order, initial value differential equations. For SU(2), a 
single complex parameter, z(t), that multiplies a + (wc 
define a± = u x ±ia y , differing by a factor of 2 from some 
others) can be solved through a single Riccati equation 
for a vector on the Bloch sphere. This z then provides 
through straightforward integration the remaining phase 
along with its breakdown into geometrical and dynami- 
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cal parts. For a spin- 1/2, charged particle in a magnetic 
field, with H = —jff ■ B(t), solution of the Schrodingcr 
equation, ill(t) = H (t)U (t) , U (0) = X, can be written 
as (an overdot will denote differentiation with respect to 
time) 



U(t) 



e z(t)a+/2 e w'(t)<T-/2 e -in(t)cr z /2 _ 



(1) 



The complex quantities z, w, and fj, are classical functions 
of time, vanishing at t = 0. Solving the Schrodingcr 
equation with this form amounts to solving a Riccati 
equation for z [1]. Unitarity of U requires: w* — 
-z*/(l + \z\ 2 ), e Im » = (1 + \z\ 2 ). Thus, there are only 
three linearly independent quantities, the real and imag- 
inary parts of z and Re /i, the last being determined by 
quadrature in terms of z. U depends on three parameters 
but the density matrix on z alone. 

A specific focus on the phase \i follows upon viewing 
the evolution operator U in Eq. (1) as U1U2, where U 2 
is the last factor involving \i and is diagonal, whereas 
U\ is given by the first two factors and depends on z 
alone. With such a product form, the Schrodingcr equa- 
tion, iU = HU, reduces to one for U 2 alone with an 
effective Hamiltonian, 



iU 2 = H cif U 2 , H eS = U^HUx - iU^Ul (2) 

This gives the equation for /x. While the factorization 
of U into two parts is generally valid, the SU(2) example 
displays explicitly the z and S 2 Bloch sphere part and the 
fi phase, respectively. Note that U\ and U 2 are not indi- 
vidually unitary, and we use the symbol tilde to signify 
this. However, as will be relevant for later development, 
factoring U 2 into real and imaginary parts, and incorpo- 
rating the former into U\ as a right multiplier gives a 
decomposition of the full U into two unitary factors, U\ 
and U 2 , the first dependent on z alone and the second 
on Re \x. These are, respectively, operations on the base 
manifold and fiber of SU(2) regarded as a fiber bundle 
[SU(2)/U(1)] x U(l). 

By defining a vector m on the Bloch sphere, 



m + 
m 3 



mi +im 2 = -2z*/(l + \z\ 2 ), 



(1 



2 )/(l + 



(3) 



with unit length: m\+m 2 + m\ = m + m- + to| = 1, one 
verifies the Bloch equation of motion, m = —2B x to, 
with a Coriolis-like appearance. The two terms HU\ 

and —iU^ 1 U\ of the effective Hamiltonian in Eq. (2) have 
equal and opposite a x and a y terms which thereby cancel. 
The effective Hamiltonian reduces to a one-dimensional 
one in a z for the fiber, the two parts, dynamical and geo- 
metrical, adding to give the full phase /x. The dynamical 
part has the expected structure of the coupling energy of 
a magnetic moment to the magnetic field, — \ih-B. At 



the end of this Letter, we cast a non-trivial problem of 
two spins into an analogous five-dimensional unit vector 
moving on the sphere S 4 according to a generalized linear 
Bloch equation. 

The above analysis extends to more general unitary Lie 
groups U, which may not necessarily be U(7V) or SU(A). 
For a sub-group B oiU, B CU, the factorization as be- 
fore, U(t) = U x {t) U 2 (t), U 2 eB, U l7 U eU, holds, along 
with Eq. (2). The factors are not unique, since for any 
b e B, we can rewrite as U — {U\b)U 2 with U 2 = b~ 1 U 2 
again an element of B. Thus, U\ is determined up to 
right-multiplication by an element b, U\b representing 
one element in the space of left cosets, U/B. Because 
b is a single phase when the group B is U(l), this non- 
uniqueness is like a gauge degree of freedom. Later below, 
we will encounter & as a full matrix operator. When U is 
SU(iV), and sub-group B is SU(A^ — 1), the above serves 
to define the fiber bundle [SU(A0/(SU(N - 1) * U(l))] 
x (SU(A r — 1) x U(l)), which is the n = 1 case of a more 
general construction [SU(AT)/(SU(iV - n) x SU(n))] x 
(SU(Ar - n) x SU(n)) with n < N/2. The base mani- 
fold within square brackets is the Grassmannian manifold 
[15] Gr c (N,n), elements of which can be represented by 
projection operators onto the n-dimcnsional subspace of 

Thus, consider the A^-dimcnsional Hamiltonian H^: 



H W = 



vt 



V 

H(«) 



(4) 



The diagonal blocks are square matrices while the off- 
diagonal V is (N — n) x n and is n x (N — n). These 
latter are taken as Hermitian adjoints and as trace- 

less for purposes of this Letter although much of our con- 
struction applies more generally. We write, analogously 
to Eq. (1), U^ft) = E/iE/b, with 



U 2 



z(t) 

l(n) 



l(N-n) 

ot 

U(JV-»)(t) 
ot U( n >(t) 



l(N-n) 

wt(t) 





I(") 



(5) 



where z and are rectangular matrices of complex pa- 
rameters. Note that Ui is constructed from block-matrix 
generalizations of a± in Eq. (1) and has blocks of zero in 
the lower and upper off-diagonal blocks of its matrix fac- 
tors. Further, U^' as a whole is unitary which leads to 
matrix renderings of the relations noted earlier between 
z and w: 



z = -7-lW, 

7l = u("-n)u(JV-»)t = I^-^+zzt, 

(7 2 ) _1 =U (n) U(") t = (iW +z t z)- 1 . (6) 

The effective Hamiltonian in Eq. (2) constructed from 
Ui is block diagonal, the off-diagonal blocks of the con- 
struction defining the equation satisfied by z, 
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HCAlV-zfV'z + HW), (7) 7 =l + z t z. (12) 



now a matrix Riccati equation [16], and the diagonal The commutator in Eq. (11) vanishes and the rest of this 
blocks giving expression is easily evaluated to give 



H, 



eff 



H(JV-«) _ z yt 

t H(")+Vtz 



(8) 



the effective Hamiltonian for the J7 2 problem. The indi- 
vidual blocks are neither Hermitian nor traceless. 

To convert Eq. (5) into a unitary decomposition, that 
is, to make both factors U\ and U2 unitary, we construct 



7i 



7 2 



(9) 



and choose its inverse square root as the "gauge factor" b 
that provides the requisite unitarization: U\~ U\b. Cor- 
respondingly, U2 = b~ 1 U2 is unitary as well. Also, using 
this U\ in an equation analogous to Eq. (2) but without 
tildes, we get the explicitly Hermitian H e g counterpart 
of Eq. (8). The upper diagonal block in such an H e g is 

^(7 1 - 1/2 ),7l /2 ]+^{7i- 1/2 (H^-")-zVt) 7 l/ 2 + h.c. } , 

(10) 

with h.c. the Hermitian conjugate of the previous expres- 
sion. The lower diagonal block of H c s is 



^|(7 2 - 1/2 ),7f] + i{7 2 - 1/2 (H(«)+ztv) 7 f + h.c.}. 

" ' " , ^ (11) 

These serve as the explicitly Hermitian effective Hamil- 
tonians for the SU(7V — n) and SU(n) problems. 

Together with the 2n(N — n) parameters in z that are 
obtained through solutions of Eq. (7), the (TV 2 — 1) pa- 
rameters of SU(iV) are thereby expressed in terms of the 
corresponding parameters in the smaller groups and one 
additional phase parameter between those two spaces 
of N — n and n dimensions. Our solution may be re- 
garded as extending the Schwinger scheme [17] which 
constructs higher spin-j representations from those of 
spin-1/2. While that scheme, and the association of 
(N = 2j + l)-level systems with spin j is familiar, it 
deals with larger representations but of the same group 
SU(2) or SO (3). Our decomposition of the full SU(7V) 
into SU(N-n) and SU(n) through two nilpotent and one 
diagonal factor in Eq. (5) exactly analogous to the sim- 
ilar Eq. (1) for SU(2) extends the scheme to a complete 
solution of an arbitrary SU(7V) Hamiltonian in Eq. (4). 

While the above constructs the fiber bundle for arbi- 
trary n, we will now concentrate on the case n = 1. The 
handling of the square root operators for the general case 
will be considered elsewhere but for n = 1, when V and 
z are column vectors, 72 in Eq. (6) reduces to a number, 



iU2,NN = [H NN + -(Vt z + riV)]U 2 ,NN, (13) 

where the lxl lower corner of the evolution operator's U2 
matrix (that is, U^ 1 '^) in Eq. (5) without the tilde) is de- 
noted by U2,nn and likewise the corner element of 
in Eq. (4) by Hnn- The equation is readily integrated 
and gives a pure phase. It complements the real factor 
(1/^/7) in the corresponding corner of Ui, these being 
analogous to the real and imaginary factors of \x for the 
SU(2) example. The geometrical contribution, —iU^U\, 
to this W-phase gives in the square bracket in Eq. (13), 
7- 1 [zt(H( Ar - 1 )-i7 JVA r)z + (ztV + Vtz)(l-7/2)]. A use- 
ful identity in these constructions is 7 = ^(V^z — z^V). 

In the effective Hamiltonian in Eq. (10) for the remain- 
ing SU(A — 1), the square roots simplify in terms of 7, 



7l /2 = I< Ar - 1 W/U/7 + l), 
7r V2 = I (Ar - 1) -zzt/(V7" + 7), (14) 



and we get 



h (at-i) = 



ZVT 



Vzt z(ztV + Vtz)zt 



V7+1 



2(y^+l) 2 



(15) 

Note that the term in square brackets in Eq. (13) is, as 
expected, — TrH^ -1 ), as a whole being traceless 

and trace being preserved in our construction. H^ -1 ) 
in Eq. (15) is the starting point for the subsequent con- 
struction of U^ -1 ) after subtracting its trace. Thus 
both the base manifold SU(iV)/(SU(iV - 1) x U(l)), de- 
termined entirely by z, and the fiber are constructed in 
a form that allows for iteration. A fiber bundle descrip- 
tion can be developed in turn for this SU(iV — 1), thus 
giving a hierarchical construction. Each step provides a 
U(l) W-phase. Turning to a geometrical picture, in the 
case of N = 2,n = 1, the Bloch sphere S 2 described 
by m in Eq. (3) is an inverse stereographic projection of 
the one-dimensional Ricmann plane given by the complex 
number z. Similarly, now for general N, and n = 1, an in- 
verse stereographic projection of the (N — l)-dimensional 
complex Riemann plane of z in Eq. (7) gives the higher 
dimensional generalization of the Bloch sphere. The vec- 
tor equation may not always reduce to a Bloch-like linear 
equation for a unit vector on S N but we now turn to an 
example for N = 4, n = 2 when it does. 

Consider two spins or qubits with a Hamiltonian of 
Spin(5) or SO (5) symmetry. There are many such within 
the full SU(4) dynamical symmetry of two spins but a 
specific representation as a concrete example is H(t) = 



4 



F 2 1(T 



(2) 

z 



(2) 



P ,r (2) 



P ^r^rr^ 



P J^rr^ 



F54«Ty , where the ten arbitrarily time-dependent coef- 
ficients F^ u (t) form a 5 x 5 antisymmetric real matrix. 
(We will use ji, v = 1 — 5 and i,j,k= 1 — 3 and summa- 
tion over repeated indices.) This form describes several 
quantum optics and multiphoton problems of four levels 
driven by time-dependent electric fields. It has more gen- 
eral couplings between the four levels than has been con- 
sidered extensively in coherent population transfer and 
other phenomena in a variety of molecular and solid state 
systems [18]. Only numerical solutions have so far been 
available but we can now provide a complete analytical 
solution. With n = 2, such a Hamiltonian in Eq. (4) has 
2x2 diagonal blocks, H^ 1 ' 2 ) = (T^fc — ^CijkFij)&k, an d 
off-diagonal V = iF^l^ + F^Oi. The matrix Riccati 
equation in Eq. (7) can be solved in terms of four real 
z M = Z4,,Zi, with z = z 4 l( 2 ) — izidi, obeying 



Zn = F 5fl (l - z v z v ) + 2F lllJ z lJ + 2F 5v z v z ll . (16) 

V and z can also be rendered in terms of quaternions 
(1,— icfi)- 7i and ~y 2 m Eq. (6) become equal and pro- 



portional to a unit matrix, (1 + z^z^)!^, and the eval- 
uation of Eq. (10) and Eq. (11) gives straightforwardly 

H^ 2) = H (1,2) - e ljk ZiF 5j a k =p F 5j Z4,aj ± F bi Ziai. Re- 
markably, in terms of a five-dimensional unit vector to, 



m ^ = (IT^' m5 =(IT^' M ' i/ = 1 - 4 ' (17) 



the nonlinear Eq. (16) in z becomes a simple, linear 
Bloch-like = 2F llv m v , fi, v = 1 — 5. As in the single 
spin case, this represents an inverse stereographic projec- 
tion, now from the four-dimensional plane z E R 4 to the 
four-sphere S* 4 , and thus provides a higher-dimensional 
polarization vector for describing such two spin problems. 
Further details and application to specific problems in 
quantum information and nuclear magnetic resonance, 
as well as other n = 1 and n = 2 applications of our 
construction, will be given elsewhere. 
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